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Abstract 

We present a general method for constructing pure-product-state representations 
for density operators of N quantum bits. If such a representation has nonnega- 
tive expansion coefficients, it provides an exphcit separable ensemble for the density 
operator. We derive the condition for separability of a mixture of the Greenberger- 
Horne-Zeilinger state with the maximally mixed state. 



1 Introduction 

The state of a quantum system composed of N subsystems is separable if it can be writ- 
ten as a classical mixture or classical ensemble of tensor product states. Although it is 
straightforward to decide whether a pure state is separable, the same question is in general 
unsolved for mixed states. A simple separability criterion is known only for two subsystems, 
and even then only if one of the subsystems is a two-level system (qubit) and the other one 
is at most a three-level system (qutrit) [0, 0] 

Recently it has been shown that there exists a neighborhood of the maximally mixed 
state in which all states are separable ||^. For qubits the size of this neighborhood 
decreases exponentially with the number of qubits 0] . These results are interesting, among 
other reasons, because NMR quantum computers |^, 0] operate with states near the 
maximally mixed state. 

In this paper we develop a general method for constructing pure-product-state represen- 
tations for A^-qubit density operators; when the expansion coefficients in such a represen- 
tation are nonnegative, it provides an explicit product ensemble for the density operator, 
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and the state is separable. In Sec. we introduce a superoperator formalism that al- 
lows us to treat different discrete and continuous representations on the same footing. In 
Sec. Ol we apply the formalism to a discrete overcomplete operator basis, leading to a 



simple rederivation of the separability bound found in Section |2.3| introduces a "canon- 
ical" continuous representation for an arbitrary density operator, and Sec. B]^ develops a 
method for deriving discrete representations from the continuous representation. In Sec. |^, 
we show that the canonical continuous representation is the unique representation contain- 
ing only scalar and vector spherical-harmonic components; all other representations can 
be obtained by adding arbitrary higher-order spherical harmonics. Section |3.2| applies the 
canonical continuous representation to mixtures of the iV-qubit "cat state" with the max- 
imally mixed state. In Sec. we derive the conditions for separability of a mixture of the 
Greenberger-Horne-Zeilinger state |^ with the maximally mixed state. Section |^ concludes 
with a brief discussion. 



2 Product-state representations and separable ensem- 
bles 

2.1 Superoperator formalism 

The set of linear operators acting on a D-dimensional Hilbert space 7i is a D^-dimensional 
complex vector space CiTi). Let us introduce operator "kets" | A) = A and "bras" {A\ = A\ 
distinguished from vector kets and bras by the use of round brackets 0. Then the natural 
inner product on CiTi), the trace-norm inner product, can be written as (v4|i?) = ti{A'^B). 
The notation S = \A){B\ defines a superoperator S acting like 

S\X) = \A){B\X) = tT{B^X)A . (1) 

Now let the set {|iVj)} constitute a (complete or overcomplete) operator basis; i.e., let 
the operator kets \Nj) span the vector space C(H). It follows that the superoperator Q 
defined by 

g^j:\Njm\ (2) 



Q,^g-'\N,) (3) 



is invertible. The operators 
form a dual basis, which gives rise to the following resolutions of the superoperator identity: 

^ = T.wm = T.\Nj)m- (4) 

An arbitrary operator A can be expanded as 

^ = E miQM) = T.NMQ]A) (5) 
j j 

and 

^ = E miNM) = j:qmn^a) . (6) 

These expansions are unique if and only if the operators Nj are linearly independent. 
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2.2 A discrete representation for the 6 cardinal directions 

The Hermitian operators Pj^, j = 1, 2, 3, /i = 0, 1, defined by 

^.. = ^(i + (-ir^.), (7) 

are the pure-state projectors corresponding to the 6 cardinal directions on the Bloch sphere. 
They form an overcomplete basis in the space of operators acting on a qubit. The super- 
operator 

^ = E \Pj,)iP^,\ = I f3|i)(i| + E (8) 

has normahzed eigenoperators \l)/^/2 and \aj)/\^, with eigenvalues 3 and 1, respectively, 
from which it follows that 

^-' = ^(^|i)(i| + EK)(^.l) • (9) 

The last equation allows us to find the dual-basis operators, 

Q,, = G-'\Pj,) = l{l + 3{-ira,) . (10) 
Any density operator p = ^{1 + S ■ a) for one qubit can be represented in the form 

p = E \Pm)iQm\p) = J E + ^i-^rs,) . (11) 

For an iV-qubit density operator, the analogous representation is 

p= w^""(ji, Pi,... ,JN, PN)Pnf.i ® ■ ■ ■ ® Pj^f^j, , (12) 

where 

w'"'''ijl,IJ-l,...jN,PN) = ir{pQj^f,^ (g) ■ ■ ■ (g)Qjj^^j^) 

= ^^'{p{^ + 3(-l)'^^or,0 ® ■ ■ • ® (l + 3(-l)'^-a,,)) . 

(13) 

This representation exists for arbitrary density operators, entangled or separable. If the 
coefficients w'^^^{ji, pi, . . . , jjy, piy) are nonnegative, the density operator p is separable, 
and the representation (|T2D provides an explicit product-state ensemble for the density 
operator. The expansions ([TT|) and (|l^) are not unique because the 6 operators P,-^ are 
not linearly independent. Throughout the paper we use the superscript "can," standing 
for canonical, to denote the natural, but generally not unique expansion coefficients that 
come from using the procedure of Sec. |2.1| . 
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For any A^-qubit density operator p, the coefficients (0) obey the bound 
ran / ■ ■ \ ^ / Smallest eigenvalue of\ 2^^~^ 

This follows from the fact that Qj^ has eigenvalues 2/3 and —1/3. The most negative 
eigenvalue of the product operator Qj^^^^ ® ■■■ ® Qjni^n is therefore (— 1/3)(2/3)^^^ = 
_22JV-i/6^, from which Eq. (|T4|) follows. 

Consider now an A^-qubit density operator that is a mixture of the maximally mixed 
density operator, 1/2^, and an arbitrary density operator pi: 



Any density operator can be written in this form 0]. The coefficients ( |T3|) for are given 

by 

which is nonnegative for 

This bound, derived in |^ using a different route, implies that, for e < 1/(1 + 2^^~^), all 
density operators of the form (|T5|) are separable. It sets a lower bound on the size of the 
separable neighborhood surrounding the maximally mixed state. 

2.3 Canonical continuous representation 

The set of all pure-state projectors on the Bloch sphere, 

P^=\n){n\ = ^{l + a-n) , (18) 
forms an overcomplete basis. The corresponding superoperator 



g = JdQ \Pn){Pri\ = vr(^|i)(i| + ^ E 



(19) 



is proportional to the superoperator of Eq. (^. The normalized eigenoperators of Q are 
1 1) / \/2 and \aj) / -\/2, with eigenvalues 2ti and 27r/3, respectively. Thus can be expressed 



as 



^-^ = ^||l)(l| + 3E|a,)(a,|| . (20) 
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The dual-basis operators are then given by 

Qn = G''\Pn) = ^il + Sa-n) . (21) 
Any density operator p = ^{1 + S ■ a) for one qubit can be represented in the form 

p = Jdn\Pn){QH\p) = ^JdnPn{l + 3S-n). (22) 
For qubits, we define the pure-product-state projector 

P{n) = Pn,(^---^Pri^ = ^(1 + a ■ ni) ® ■ ■ ■ ® (1 + a ■ nTv) , (23) 

where n stands for the collection of unit vectors ni, . . . , riN. Any A^-qubit density operator 
can be expanded as 



P 

where 



J dQn w''^-^{n)P{n) = J dQn, ■ ■ ■ dnn^w'^^{n)P{fi) , (24) 



w'='^'^(n) = tr(pg(n)) (25) 
is the canonical expansion function, with 

Q{n) = Qr,^®---®Qn^ = J^i^ + 3a ■ ni) ® ■ ■ • ® (1 + 3a ■ n^v) • (26) 

If the coefficients w'^^'^{n) are nonnegative, the density operator p is separable, and the 
representation ( pij ) provides an explicit product-state ensemble for the density operator. 

Using the same argument as in the preceding subsection, one can show that for any 
density operator p 

ran/~\ ^ / smallcst elgcuvalue \ 2^^~^ 



Applying this result to density operators of the form ([Iq), one arrives again at the lower 
bound (p!?!) on the size of the separable neighborhood of the maximally mixed density 
operator 

2.4 Other discrete representations 



The expansion coefficients (0) for the discrete representation of Sec. |2]^ can be obtained 



by evaluating the expansion coefficients ( pSf ) for the canonical continuous representation 
along the cardinal directions for each qubit and then renormalizing the resulting function. 
This is no accident. It is easy to characterize a class of discrete representations whose 
expansion coefficients are similarly related to the continuous representation. 
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Consider a discrete set of projection operators, 



-Pn, = = ^(1 + (?• no) , a = l,...,K. (28) 



The corresponding superoperator 



K 



a=l 



\ [k\i){i\ + ■ + ■ + E El 



na)j{na)k (29) 



generates dual-basis operators and expansion coefficients proportional to those for the con- 
tinuous representation if and only if Q is proportional to the superoperator (p!9| ) for the 
continuous representation, i.e., if and only if 

= E^"' (30) 

a 

1 1 

o^Jk = ^EK)iK)fc- (31) 
6 K ^ 

The trace of Eq. (|3l|) is satisfied by any collection of unit vectors. 

When these conditions are satisfied, the superoperator ( pQ] ) simplifies to 



with an inverse 



^ = f (|i)(i| + ^El^.)(^.l) > 
^-^ = ^(|i)(i| + 3Eh)(^.l) 



(32) 



(33) 



that generates dual-basis operators 

Qn. = G-'\PnJ = ^(1 + 3a ■ n„) . (34) 
Any density operator p = i(l + S* ■ a) for one qubit can be represented in the form 

p = ^|P„)(g„|p) = ^ ^P„(l + 35^-n„) . (35) 



The conditions (|30| ) and (|31| ) are satisfied by unit vectors that point to the vertices of 
any regular polyhedron inscribed within the Bloch sphere — i.e., a tetrahedron, octahedron, 
cube, icosahedron, or dodecahedron. The 4 projectors for a tetrahedron are linearly inde- 
pendent, making the corresponding tetrahedral representation ( ^51) unique. An octahedron 
gives rise to the 6-cardinal-direction representation of Sec. |2.2| . The regular polyhedra by no 
means exhaust the possibilities for representations of this sort. One can use simultaneously 
the vertices from any number of polyhedra. Moreover, one can produce an appropriate set 
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of unit vectors by the following procedure: starting with any set of unit vectors in the first 
octant, reflect the set successively through the three Cartesian planes to produce a set of 
vectors that occupies all eight octants. 

For more than one qubit, the possibilities for discrete representations that follow from 
the continuous representation are even more diverse. The simplest possibility is to use the 
same set of pure-state projectors for each qubit, but this is not necessary. One can use a 
different set of projectors for each qubit, or more generally, one can use a "tree structure," 
in which there is a different set of projectors for the nth qubit for each projector for 
the first n — 1 qubits. There are possibilities more general even than those with a tree 



structure, analogous to the "domino states" introduced for two qutrits in |T0[. We do not 
dwell here on the myriad of discrete representations of this sort, merely noting that since 
the expansion coefficients in such representations follow from the expansion coefficients 
for the continuous representation, all such representations lead to the same general lower 
bound (|T^ on the size of the separable neighborhood surrounding the maximally mixed 
state. Both the continuous and these discrete representations can, however, give better 
bounds on the separability of particular states of the form (^), a point to which we return 



in Sec. 3.2 



3 Spherical-harmonic expansions 

In this section we explore the relationship between product-state representations and sphe- 
rical-harmonic expansions. 

3.1 Spherical-harmonic expansions and the PauU representation 

We begin by noting that if is a unit vector, then 

a ■ n = CT3 cos 6 + ai sin 6 cos + (72 sin 6 sin 

^(^a^r + -i^^i-Yr + Y.-n - -2^(n^^* + Y,-n) (36) 



Yo=^, Y,' = Jocose, r±i = W^sin^e±'^ (37) 



where 



are spherical harmonics and, in the last line, we define 

a° = v^l, a? = ^(73, a^' = T^icj^±ta2). (38) 
We can then write any pure-state projector on the Bloch sphere in the form 

Pn =\n){n\ = l{l + a.n) = ^j:j: aTYr • (39) 

;=o m=—l 
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Any A^-qubit density operator can be written in the form 

p= f dVLnw{fi)P{n) , (40) 



where the expansion coefficients w{n) = w{ni, . . . , f?jv) are not necessarily nonnegative [the 
expansion coefficients w(n) need not be those of the canonical representation (|25|)]. A 
density operator p is separable if and only if it can be written in the form ( ^OD for some 
nonnegative probability density w{h). Now expand the product projector (^) in terms of 
spherical harmonics, 

pin) = ^ E >^r--->^r*<^ (41) 

{ij<l,mj} 

where 

1 h 1 In 

E - E E ■■■E E (42) 

{Zj<l,r?ij} li=0 mi=~li (jv=0 m]v=— ijv 

denotes a sum that contains only scalar and vector terms (/ = 0, 1), and also expand the 
expansion coefficients, 

M^)= E «Dr^r--->^r> (43) 

where 

00/1 CO In 

E - E E ■■■E E (44) 

{Ij.rrij} li=0 mi=—li i]v=0 m]v=— Zjv 

denotes an unrestricted sum over spherical harmonics. Inserting Eqs. (^Tj) and ( ^3]) into 



Eq. (|40D and using the orthonormality of the spherical harmonics yields the (unique) Pauli- 



operator representation of p: 



P=^ E (45) 

{Zj <l,mj } 



We see that p depends on and determines only the scalar and vector parts of a^^ 



mi ...mjv . 

■ In ' 



the higher-order-spherical-harmonic (HOSH) content of w{n) corresponds to the freedom 
in constructing pure-product-state representations of p. It is straightforward to generate 



all pure-product-state representations: given the Pauli-operator representation (|45|) of p, 
any w{n) is given by Eq. (^), where the coefficients ci]^^"i^^ are arbitrary except for 



/i, . . . , /at < 1. 

The canonical continuous representation w'^'^^{n) of Eq. (^5)) is the unique pure-product- 
state representation given by the / = and / = 1 components. We can see this directly in 
the following way. Write the Pauli representation of p in the form 

P= E c„,...„^a„, ® • • ■ ® c^ajv . (46) 

ai,...,ajv 

where 

Cai...aN = tr(pcra^ ® ■ ■ ■ (g) (T„j^) = ((Ta, (g) • ■ ■ (g) a^iy) ■ (47) 
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Here cxo = 1, and the sums over Greek indices run from to 3. Now we have that 



w^^''{n) = -^^^((l + Sa-ni)®---® (l + 3CT-n7v; 
3 \^ 

7-) ((^l)aiO-ai ® • ■ ■ ® (^Ar)ajvC^ajv) (48) 

3 \^ 



where {nj)a. = 1/3 if = and {nj)aj is a Cartesian component of if = 1,2, or 
3. Thus any pure-product-state representation of a density operator p can be obtained 
by adding HOSH to the canonical continuous representation i(7^^"(fi). The discrete repre- 
sentations discussed in Sec. are examples — though by no means the only examples — of 
representations with HOSH content. 

The continuous canonical representation can be generated directly from the Pauli rep- 
resentation of p by making the substitutions 

2^p — y w=^'^(n) , (49) 
These rules are equivalent to cr™ — > F^"^ and to cTq, — > (3/47r)nQ,. 



3.2 The e-cat state 

As an example, consider the following A^-qubit generalization of the Werner [|ll| state, 

Pe = + e|^cat)(^/'cat| , (51) 

where 

|V^cat) = ^(|0...0) + 11...1)) (52) 

is the A^-qubit "cat state." We call the mixed state (|51| ) the e-cat state. 

Now let |0)and |1) be the eigenstates of (T3, with |0) at the top of the Bloch sphere and |1) 
at the bottom, so that ai = |0) (1| + |1) (0|, da = -i|0)(l| +2|1)(0|, and = |0)(0| - |1)(1|. 
It follows that the Pauli representation of the e-cat state is 

= ^((1-^)1 

+ |(1 + dg) ® ■ ■ ■ ® (1 + (T3) + |(1 - dg) ® ■ ■ ■ ® (1 - dg) (53) 



+ -(CTl + 10-2 



) ■ ■ ■ {ai + icTa) + ^(cTi - icra) ® ■ ■ ■ ® (di - 2(12)^ . 
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Applying the rules (|i9| ) and (0) yields the continuous canonical representation of p^: 

+ ^(1 + 3 cos^i) • ■ ■ (1 + 3 cos^jv) + ^(1 - 3 cos^i) ■ ■ ■ (1 - 3 cos^jv) 

+ e 3^ sin 6^1 ■ ■ ■ sin On cos(0i + . . . + 0iv)^ . (54) 

The minimum value of wl^'^{n) occurs either (i) when one qubit is evaluated at the 
north (south) pole, and the other N — 1 qubits are evaluated at the south (north) pole, or 
(ii) when all the qubits are evaluated on the equator with cos(0i + . . . + ^^v) = —1. For 
N = 2 and = 4, both the poles and the equator give the same minimum value, which 
shows that the e-cat state is separable for e< l/9foriV = 2 and for e < 1/81 for = 4. 
For A^ = 3 and A^ = 5, the equator gives the minimum value, which shows that the e-cat 
state is separable for e < 1/27 for A^ = 3 and for e < 1/243 for A^ = 5. For A^ > 6, the 
poles give the minimum value, which shows that is separable for 

1 

^ — 1 ± 2^ + 22^-2 ' ^ ^ 

the upper (lower) sign applying when A^ is even (odd). Notice that except for N = 2, all 
these bounds are better than the general-purpose bound of Eq. (|T7|). 



These bounds on the separability of the e-cat state can be improved by using tetrahedral 



ensembles of the sort discussed in Sec. |2.4| . The tetrahedral ensembles are chosen to have a 
tree structure such that the vertices avoid the places where wf'^ln) has its smallest values. 
We do not discuss these tetrahedral bounds here, preferring instead to consider in detail 
the separability of the e-GHZ state. 



4 The e-GHZ state 



In this section we consider the 3-qubit version of the e-cat state introduced in Sec. |372 . 
In the case of 3 qubits, the cat state of Eq. (0) is called the Greenberger-Horne-Zeilinger 
(GHZ) state [|, 

|V^GHz) = ^(|000) + |lll)) ; (56) 

we call the mixture of Eq. (|5TD the e-GHZ state. We show in this section that the e-GHZ 
state is separable if and only if e < 1/5. The class of e-GHZ states has been considered 
previously in [0. 

4.1 Separability of the Werner state 

Before proceeding to the proof of the separability of the e-GHZ state, we review a similar, 
but simpler proof for the Werner state [|1T|, the 2-qubit version of the e-cat state: 

Pw = ^1®1 + |(|00) + |11))((00| + (11 
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2 (l ® 1 + e (0-3 ® 0-3 + (Ti (g) CTi - (72 ® (T2)) . (57) 



The Werner state is separable if and only if e < 1/3 ||T3|, |T3|, a result we now show in the 
following way. 

A density operator p for 2 qubits, A and B, is separable if and only if it can be written 

as 



P 



with w{nA,nB) being a nonnegative probability distribution. If a state is separable, the 
coefficients (^7|) in the Pauli representation (|46|) of p become classical expectations over the 
distribution w{nA,nB), expectations that we denote by E[-]: 



Caf3 = tr(p aa = j dVLfi^dVt^j^ {nA)a{nB)(} w{nA, Ub) = E[{nA)ainB)f3] ■ (59) 

For the Werner state, the coefficients Cq,/? can be read off the Pauli representation in 
Eq. (0): 

E[inA),] = E[{nB)j] = , (60) 
E[{nA)jinB)k] = , hij^k, (61) 
E[{nA)iinB)i] = -E[{nA)2inB)2] = ^[(n^)3(riB)3] = e • (62) 

If w{nA,nB) is a probability distribution, we have that 

E[{nA)j{nB),]\ < ^ElinAY, + (n^)?] • (63) 
Adding over the three values of j gives 

3e < \E[inA)l + Ml + Ml + (n^)? + Ml + Ml] = 1 , (64) 

which shows that if e > 1/3, the Werner state is nonseparable. 

We still have the task of designing a product ensemble that gives the Werner state for 
e < 1/3. Defining three density operators 

Pi = ^{Pe, ® Pes + P-e, ® P-es) = ^{l ® I + ^3 a^) , (65) 

1 1 

P2 = ^{Pe, ® Pe^ + P-e^ ® P-e^) = -(1 ® 1 + fTl ® fTl) , (66) 

P3 = ^{Pg, ^ P-e-2 + P-e, ^ Pe,) = ^{l ® I - ^ , (67) 

we see that each takes care of one of the correlations in Eq. (|62l). An equal mixture gives 
the e = 1/3 Werner state: 

1 1/1 \ 

^(Pl + P2 + Ps) = ^ ( 1 ® 1 + ^{Cr3 O 0-3 + (Ti (g) (Ti - 0-2 (S) ^2) j . (68) 

Smaller values of e can be handled by adding on a product ensemble for the maximally 
mixed state. 
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4.2 Separability of the e-GHZ state 

Consider now the e-GHZ state 
1 -e 



PeGHZ 



1® 1 + -(|000) + |111))((000| + (111|) 
. (g) 1 (g) 1 + e(^l (g) 0-3 (g) 0-3 + 0-3 ® 1 (g) 0-3 + 0-3 (g) 0-3 ® 1 

+ (Ji (g (Ji (g (Ji — (Ji (g (72 ® 0-2 — 0-2 ® O"! ® — ® 0-2 ® CTl) 



(69) 



A density operator p for 3 qubits, A, 5, and C, is separable if and only if it can be 
written as 

p = j dVtfij,dVtfigdVtn^ w{nA, Ub, nc)PnA ® P-Hb ® ^nc ' C^O) 

with w{nA,nB,nc) being a nonnegative probability distribution. If a state is separable, 
the coefficients (^) in the Pauli representation of p become classical expectations over 
the distribution w{nA,nB,nc)- 

= tr(paa (g (g (J^) = E[{nA)a{nB) f3{nc)-y] . (71) 

For the e-GHZ state, the coefficients Ca/sj can be read off the Pauli representation in 

Eq. my- 

E[{nA)j] = E[{nB),] = E[{nc)j] = , (72) 

E[{nA)j{nB)k] = E[{nA)j{nc)k\ = E[{nB)j{nc)k] = e^js^s , (73) 

E[{nA)j{nB)k{nc)i] = e (^ji^i^n - 5ji425/2 - 5j24i5«2 - 5i242^«i) • (74) 

Now define a "vector" 

N = [{nB)i{nc)i - (ris) 2(^0)2)61 - ((^5)1(720)2 + (nij)2(nc)i)e2 + (^1^)363 , (75) 

whose magnitude satisfies 

N-N = ((nB)i(nc)i - (725)2(^0)2) + {inB)iinc)2 + inB)2{nc)i) +(^5)3 

= {Ml + Ml) {inc)l + inc)i} + Ml (76) 
< UB-nB = 1 ■ 

If w{nA, ubi nc) is a probability distribution, we have that 



5e 



E 



UA-N 



< E 



UA-N 



< -E 
- 2 



ua- nA + N ■ N 



< 1 



(77) 



which shows that if e > 1/5, the e-GHZ state is nonseparable. 
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To finish the proof, we need to construct a product ensemble for the e-GHZ state when 
e < 1/5. We begin by defining five density operators 

= ^(l® l®l + l®cr3®cT3 + (j3(g)l®a3 + cr3®cr3®l) , (78) 

P2 = \ {Pe, ® Pe, ® Pe, + Pe, ® P-e^ ® P~e^ + P~e^ ® Pe^ ® P-e, + P-e^ ® P-e^ ® Pe^ 

= ^(l®l®l + ai(8)(Ti®(Ti) , (79) 
8 ^ 

P3 = ^ {Pe, ® Pe2 ® P-e2 + Pei ® P~e2 ® + i^-el ® Pe^ ® ^e^ + P-el ® P-e^ ® P-e^^ 

= ^(1(8)1® 1 -ai (g) a2®f72) , (80) 
8 ^ 



^(l®l®l-a2®(Ti®(T2) 



^11 



= ^(l® I®l-a2®a2®ai) , (82) 

each of which takes care of one of the nonzero correlations in Eqs. (|73D and (|74D. An equal 
mixture gives the (e = l/5)-GHZ state: 

^(Pi +P2 + P3 + P4 + P5) 
5 

= ^ f 1 ® 1 ® 1 + 1 ® Cr3 ® ^3 + ^3 ® 1 ® (T3 + (T3 ® (73 ® 1 (83) 

8 V 5 ^ 

+ ai ® (Ji ® (Ji - (Ji (g) a2 (8) 0-2 - 0-2 ® 0-1 ® 0-2 - 0-2 (8) 0-2 (g) (Ji) ^ , 

thus completing the proof. 

For both the Werner state and the e-GHZ states, we are able to construct an opti- 
mal ensemble based on the 6 cardinal directions, but notice that this ensemble is not the 



canonical representation of Sec. 



5 Conclusion 

In this paper we introduce a canonical procedure for constructing representations of den- 
sity operators in terms of complete or overcomplete sets of pure product projectors. The 
procedure is useful for investigating the separability of density operators: if the expansion 
coefficients in such a representation are nonnegative, they can be interpreted as a probabil- 
ity distribution, with the result that the representation then provides an explicit product 
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ensemble for the density operator. Generally, however, the canonical representations have 
negative values even when the density operator is separable. 

The canonical continuous representation is the unique representation that contains only 
scalar and vector spherical-harmonic terms, and we outline how to obtain all possible pure- 
product-state representations from the canonical continuous representation by adding arbi- 
trary higher-order-spherical-harmonic content. Searching the space of spherical-harmonic 
expansions for nonnegative representations would be quite difficult, but we provide an opti- 
mal ensemble in the particular case of the e-GHZ state, thereby identifying the separability 
boundary for this class of states. 
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Note added 

The separability bound for the e-GHZ state derived in Sec. ^ has also been found by Diir, 
Cirac, and Tarrach |T5[ as part of a more general result. They prove that the A^-qubit e-cat 



state (pID is separable if and only if e < 1/(1 -|- 2 ), which supersedes our Eq, 
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